Spin Exciton Formation inside the Hidden Order Phase of CeBg 
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The heavy fermion metal CeBe exhibits hidden order of antiferroquadrupolar (AFQ) type below 
Tq — 3.2K and subsequent antiferromagnetic (AFM) order at Tn — 2.3K. It was interpreted as 
ordering of the quadrupole and dipole moments of a Fs quartet of localised Ce 4/^ electrons. This 
established picture has been profoundly shaken by recent inelastic neutron scattering [l] that found 
the evolution of a feedback spin exciton resonance within the hidden order phase at the AFQ wave 
vector which is stabilized by the AFM order. We develop an alternative theory based on a fourfold 
degenerate Anderson lattice model, including both order parameters as particle-hole condensates of 
itinerant heavy quasiparticles. This explains in a natural way the appearance of the spin exciton 
resonance and the momentum dependence of its spectral weight, in particular around the AFQ 
vector and its rapid disappearance in the disordered phase. Analogies to the feedback effect in 
unconventional heavy fermion superconductors are pointed out. 



PACS numbers: 71.27.+a, 75.30.Mb, 75.40.Gb 

In strongly correlated f-electron metals the investi- 
gation of hidden order (HO) of unconventional non- 
magnetic type is a topic of central importance [2 . The 
most prominent and most investigated heavy fermion 
compounds that exhibit HO at low temperatures are 
URu2Si2 and CeBg which have tetragonal (D4/J,) or cubic 
{Oh) structure respectively. Two issues arise in the con- 
text of hidden order: Firstly, which symmetry is broken 
in the HO phase and to which irreducible representation 
the order parameter belongs. Secondly, should the or- 
dering be described as appearance of spontaneous long 
range correlation between local f-electron degrees of free- 
doms, i.e. f-electron multipoles, or should HO rather be 
described as condensation of itinerant heavy particle-hole 
pairs with a nontrivial orbital structure. These opposite 
perspectives have prevented a clear identification of the 
HO in URu2Si2 until present. 

On the other hand since the work of Ohkawa[3] the HO 
in CeBg which appears at Tq=3.2 K has always been 
taken granted as a paradigm of the localised HO picture. 
In subsequent work along this line[4l |5] it was clarified 
that the primary HO parameter is of the two-sublattice 
antiferroquadrupolar (AFQ) F^ type {Oyz^Ozx^Oxy) 
with wave vector = (^, ^, ^) in r.l.u.(R-point) which is 
nearly degenerate with an antiferrooctupolar (AFO) F^ 
i^^xyz) order parameter which is strongly induced in an 
external field. Here ± denotes the parity with respect to 
time reversal. The hidden multipolar order parameters 
are supported by the fourfold degenerate 4f crystalline 
electric field (CEF) ground state Fg. This localized sce- 
nario explains a large body of experimental results, in- 
cluding the field dependent increase and anisotropy of the 
critical temperature and the field induced Bragg peaks [6] 
at and NMR line shifts [5^ although there is no macro- 
scopic symmetry breaking observed f^. A further im- 
portant support for this picture comes from the predicted 
rapid field induced increase of the secondary octupole or- 



der parameter [4] which was directly confirmed by RXS 
experiments [8]. At temperatures below Tat =2. 3 K fi- 
nally CeBe develops antiferromagnetism (AFM) with 
Q = (i^ 0) (^ S -point) that coexists with AFQ 
order. Important information on HO may also be gained 
from the magnetic excitation spectrum. For finite fields 
that stabilizes the AFQ/ AFO HO it was investigated 
within generalized Holstein-Primakoff and random phase 
approximation (RPA) approaches [9| |10]. Both lead to 
multipolar excitation bands in the range 1 — 2.5 meV 
and for finite applied field [11 their salient features agree 
with experimental results from inelastic neutron scatter- 
ing (INS). In the numerous theoretical investigations of 
HO in CeBe the localised 4f approach was chosen and 
itinerant 4f character was completely neglected. This 
seems surprising because CeBe is a prominent exam- 
ple of a heavy fermion metal with one of the heaviest 
masses reported (m*/m > 17) [12 and the Ce- dilute La- 
substitutes [13 being the standard case of Kondo reso- 
nance dominated local Fermi liquids with all the typical 
Kondo anomalies identified there. In fact the estimated 
Kondo temperature of the concentrated CeBe from 
quasielastic neutron scattering[14 is T* c:^ 4.5 K which is 
of the same order as Tq and T^. Therefore one question 
is whether the HO physics of CeBe can be completely 
explained within the conventional localized 4f approach. 

Recent zero field high resolution INS experiments by 
Friemel et al[l have indeed seriously questioned the stan- 
dard picture and found intriguing new evidence that 
the dynamical magnetic response in the HO phase can- 
not be understood in the localized approach and, as in 
URu2Si2, requires taking into account the itinerant quasi- 
particle nature of f electrons. It was found that the low 
temperature magnetic response within the HO phase is 
determined by a pronounced feedback effect^ i.e. a modifi- 
cation of magnetic spectral properties due to the appear- 
ance of order parameters: i) Below T/v a spin gap opens 
for low energies and spectral weight from the quasielastic 



2 




FIG. 1. (Color online) Non-interacting susceptibility at the R-point (|, |, |): (a) imaginary part and (b) real part. Inset of (a) 
shows the quasiparticle DOS in PM and coexisting AFQ/AFM (T=6 x 10~^t) phase where fi — — 0.06t is the chemical potential. 
Inset of (b) gives the schematic temperature dependence of order parameters, quasiparticle model parameters: t = 22.4 meV; 
V = 0.3t; £/ = — O.Olt; gap parameters: Aq/ = 0.015t; Aq/ = 0.005t (c) Imaginary part of RPA susceptibility at R-point (inset 
shows the model for the quasiparticle interaction Jq along FR direction with Jq/ = O.lt). 



region [15 is shifted to higher energies forming a pro- 
nounced resonance at with peak position uor — 0.5 
meV. ii) Using the single-particle charge gap 2 A :^ 1.2 
meV in the HO phase from point-charge spectroscopy p!6l 
00^^/2 A = 0.42 < 1 is fulfilled showing that the resonance 
is indeed split off from the continuum, iii) The reso- 
nance appears mainly at the AFQ but not at the 
AFM Q vector and shows no dispersion. Its intensity 
decreases rapidly when approaching T^r from below in 
an order-parameter like fashion. These characteristics 
of the magnetic spectrum in CeBe do not suggest a 
spin wave origin but rather are reminiscent of spin ex- 
citon resonances observed before in Fe-pnictide[17 and 
heavy fermion superconductors [18j [19] as well as Kondo 
insulators [20*, ^1 . The results of Ref."l are the first 
clearcut example of the feedback spin exciton appearing 
within the AFQ HO phase. This proves that the local- 
ized 4f-scenario for CeBe is not adequate to explain its 
intriguing low energy spin dynamics and its momentum 
dependence. 

In this Letter we therefore propose and explore an al- 
ternative route of theoretical modeling. We start from 
the central idea that the AFQ and AFM order parame- 
ters are to be described as particle hole condensates in the 
itinerant heavy quasiparticle picture. The latter is ob- 
tained from a microscopic fourfold (Fg-type) degenerate 
Anderson lattice model. It includes both twofold pseudo- 
spin (cr =1, \) and twofold pseudo-orbital (r = ±) degen- 
eracies of the hybridizing conduction (c) and 4f electron 
(f) in the Fg CEF ground state according to 

^ = [^k^L^km + 4/km/km + (4,m/km + h.c}j 

k,m 

+ E UjffLhnflhm. (1) 

i,m,n 

where m = (r, a) represents the fourfold Fg degeneracy. 
Here cj^^ creates an conduction electron in the chan- 
nel with corresponding Fg symmetry and wave vector 
k. Furthermore, and ej^ = are effective tight bind- 
ing dispersions of the conduction band and the atomic 



/ level position respectively. For the former we restrict 
to the next neighbor hopping (t), i.e., eg. = 2t^^ cos kn 
{n = x^y^z) which leads naturally to the AFQ ordering 
vector Q^ Furthermore creates the f electron with 
momentum k, and Uff is its on-site Coulomb repulsion. 
Finally Vk is the hybridization energy between the low- 
est 4f doublet and conduction bands which contains in 
principle the effect of spin orbit and CEF but is taken as 
constant Vi^ = V here. 

In the limit Uff ^ oo double occupation of the f- 
states are excluded, this is achieved by using the auxil- 
iary boson bi at each site i, with the constraint blbi + 
Y^mfimfim = 1- ^hc mean field (MF) approxima- 
tion (r = {bi) = {b])) diagonalization leads to hybridized 
quasiparticle bands[13j. They are determined by renor- 
malized f level ^j^ = ej^ + A and effective (reduced) hy- 
bridization Vk = rVk- Minimizing the MF ground state 
energy leads to selfconsistent equations for r, A. 

The AFQ and AFM order parameters with wave vec- 
tors Q and respectively contribute extra MF terms 

y^AFQ = X]^Q'(/k,+a/k+Q'-a + /k _a/k+Q',+a), 
kcr 

y^AFM = ^ ^Qiflr^h+Qrl + /kri/k+Qrt)' (2) 
kr 

Our emphasis in this work is on the feedback effect, i.e. 
the effect of the gap opening within the HO phase on 
the magnetic response. Therefore we do not attempt a 
microscopic calculation to derive these order parameters 
and their temperature dependence. We include them as 
symmetry breaking molecular field terms in the Hamil- 
tonian and take a generic empirical temperature depen- 
dence. The MF Hamiltonian Hmf obtained from Eq. ([T]) 
is diagonalized by the unitary transformation 

/km = 'W+,k«+,km + '^-,k<^-,km 
Ckm = '^-,k<^+,km " '^+,k<^-,km- (3) 

where 2ul^ = 1 ± (e£ - et)/\/(4 " 4)' + 4^4' > ^ad- 
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ing to Hmf = E -£^k«l km«",km + -^(r^ - 1), where 



Ei = i[4+ek±V - ^t)' + 4^k] are the pair (a = ±) 
of hybridized quasiparticle (ao^km) bands, each fourfold 
(m=l-4) degenerate. Here = V^{1 — Uf) denotes the 
effective hybridization obtained by projecting out double 
occupancies. Due to 1 — n/ <C 1 Vk is strongly reduced 
with respect to the single particle Vk which leads to the 
large quasiparticle mass. Introducing new Nambu opera- 
tors as 4 = (ct, C^+q,, C^+q) where Cl = {bl^^, bl^) 

and 6^ = (al,k+^,al,k+;,al,k_^,al,k_;), we can write 
the total Hamiltonian Htot = ^mf + ^afq + ^afm as 



^tot ^^k/^kV^k; /5k 



^k Aq. Aq 
Aq/ £^k+Q' 
Aq £^k+Q, 



here £^k = ^k (8) tq ctq, Aq/ = AQ/(pk,Q' (8) fo (8) a^) 
, and Aq = AQ(pk,Q <S) <S) ^o), where fk and pk,Q/ 
are 2x2 matrices in a = ± space with matrix elements 
^k^ = ^a/3^k' and p^^, = iiak^/3,k+k'. cr/,T/ are the 
Pauli matrices acting in pseudo-spin and pseudo-orbital 
space, respectively. 

Defining the Matsubara Greens function (GF) ma- 
trix as G\^{r) = — (T?/ik(T)?/i^(0)), and solving the stan- 
dard equations of motion, one can find G\^{uJn) = 

(^iujn — /^k^ which can be written as 



Gk(c^n) 



AO 

'-^k+Q',k 
'-^k+Qk 



AO 

'-^k,k+Q' 

AO 

'-^k+Q' 

AO 

'-^k+Q,k+Q' 



AO 

'-^k,k+Q 

AO 

'-^k+Q',k+Q 

AO 

^k+Q 



,(4) 



here is a 8 x 8 Green's function matrix in (a, m) 
space. For the magnetic excitation spectrum we need 
the dipolar susceptibility matrix given by Xq (^) = 



XI /k+qm^mm'/km' 



kmm' 



-0{t){Tf^{t)£^{O)), where 

are the physical magnetic dipole operators I' = x, z). 
In cubic symmetry it is sufficient to calculate X^{^)^ 
corresponding toJW = |fo o-^, defining s = 

(a, k + q, mi) and s' = {a' ^ k, 1712) one finds 



Xo(q,^)= Xq^(^) 2^ {pI 

a.oi'\<im-\_m2 

GlAv + u:')Gl,A^')U^^+^o^ (5) 

Here the ^ are the matrix elements of reconstructed 
quasiparticle states in the AFQ/AFM state. They play 
a similar role as the 'coherence factors' in the spin exci- 
ton formation in unconventional superconductors. The 
dynamic magnetic susceptibility in RPA has the form 



XRPa{^,Oj) = [1 - JqXq"(w)]"'Xq"(w), 



(6) 



where Jq is the heavy quasiparticle interaction taken 
diagonal in (a, m) band indices. In principle it is 
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FIG. 2. (Color online) Contour plot of imaginary part of 
RPA dynamical susceptibility (a) from r(0 0) to R(| | |); 
(b) from X(0 |) to R-point; (c) from X-point to R-point; (d) 
from A(0 0|)toT(| | |) (note different scale). Resonance 
is located around R and shows little dispersion. 



determined by processes beyond the slave boson MF 
approximation[22j. However as in other spin-exciton 
theories [21] |23] we adopt here an empirical form of 
Lorentzian type that is peaked at the AFQ ordering vec- 
tor where the resonance appears. 

We will now discuss the characteristics of the magnetic 
excitation spectrum obtained from x!rpa(^'^) show 
that it explains all the essential experimental features 
observed in CeBe. In accordance with the heavy quasi- 
particle mass in this compound the chemical potential is 
chosen close to the top of the lower quasiparticle band 
(Fig. [1^ inset: /i = — 0.06t) where dispersion is ffat, lead- 
ing to a realistic mass enhancement m*/m 20. All 
other model parameters are deffned in Fig. [l] 

First the spectrum Xq{c\^uj) of non-interacting quasi- 
particles is shown in Fig{l^ with constant-q scans for 
the paramagnetic (PM), AFQ and coexistent AFQ/AFM 
phases, respectively. In the PM state the spectrum ex- 
hibits the cf- hybridization gap at the R-point. When 
the AFQ, AFM order appears their corresponding gaps 
Aq/ and Aq push the magnetic response to higher en- 
ergies. The associated real part in Figjl]3 then shows a 
much enhanced response at these energies. As a conse- 
quence the magnetic spectrum for the interacting quasi- 
particles may develop a resonance when the real part of 
the denominator in Eq. ([6| is driven to zero equivalent to 
a pole in Xrpa{Q! -t^)- Due to the 3D electronic struc- 
ture Xo(Q^ ^) will not be singular and the resonance will 
only appear for J(Q') larger than a threshold value. The 
imaginary part is generally non-zero but small leading 
to a large resonant response at the pole position. The 
resonance appears in the HO phase when Jq/ / 1 lies in a 
reasonable range such that the pole exists only when the 
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FIG. 3. (Color online) Contour plot of imaginary part of 
RPA dynamical susceptibility in (lilil)-plane of the reciprocal 
space; (a) at cj = cjr = 0.07t (spin exciton resonance energy) 
pronounced localized peak at R-point appears (b) for energy 
in the spin gap, i.e., cj = O.ScJr- Intensity at R-point vanishes 
due to spin gap formation. 



real part is enhanced by the gap formation. Then the res- 
onance condition 1 = JQ'Xo(Q^ ^r) is fulfilled only in the 
AFQ ordered regime. The magnetic spectrum of inter- 
acting quasiparticles is shown in Fig{T]3. It shows indeed 
a peak appearing in the AFQ phase and a sharp reso- 
nant peak at uJr/^Ac = 0.64 at low temperature when 
both gaps are present. Here Ac = 0.056t is the charge 
gap given in the inset of Fig{l^. This explains the central 
observation of the R-point resonance in CeBg. 

The momentum dependence of the spectrum in the 
AFQ/AFM phase and in particular the resonance peak 
is shown in Fig. [2] as contour plot in the q, uj plane with 
the wave vector q chosen along various symmetry direc- 
tions. There are two main characteristics: i) The single- 
particle spin gap due to the hybridization and enhanced 
by the AFQ/AFM gap formation appears most promi- 
nently close to the R point and less at other symme- 
try points like, e.g., T(^, ^, |). ii) The many-body res- 
onance peak is also strongly constrained to the narrow 
region around the R-point, partly due to the suppression 
of the Xo(q, ^) peak (Fig. [T]3) when q moves away from 
i?(^, ^, ^) and partly due to the decrease of Jq. Both 
mean that the above condition for the resonance can only 
be fulfilled in a narrow region around the R-point where 
it is almost dispersionless. This corresponds exactly to 
the experimental observation in CeBg and similar ob- 
servations have been made in the Ce-based superconduc- 
tors [ISl H^. A complementary constant-cj plot of the 
magnetic scattering intensity which is proportional to 
x!rpa(^'^ = const) is shown in Fig. [s] for q in the (hhl) 
plane as in the experimental scattering geometry. At the 
resonance position oor (a) the momentum dependent scat- 
tering intensity is strongly peaked at the R-point with 
rapid decay in all q directions into the scattering plane. 
On the other hand for uj = O.SuJr (b) which is in the spin 
gap region, the latter shows up as a complete depletion 
of intensity at the R-point. Due to the magnetic sum 
rule the formation of the spin gap at this energy leads 
to a roughly even redistribution of the spectral weight 
across the whole scattering plane. This complete change 
of constant-cj intensity in (hhl) plane for uj = ujr and 



<C cj^ is in agreement with the experimental result [T]. 

Now we discuss the temperature dependence of reso- 
nance intensity. We start from itinerant type AFQ/AFM 
order parameters in Eq. ^ with a typical MF BCS tem- 
perature dependence shown in Fig.[T|3 (inset). The reso- 
nance intensity at the HO wave vector in Fig. [T]^ appears 
already at Tq and is further enhanced below T^v- Ex- 
perimentally it is found that it is strongly suppressed in 
the region Tn <T < Tq. This is an effect of quadrupole 
OP fluctuations at zero field due to the near degener- 
acy with octupole order [24 which strongly suppress its 
amplitude. For example specific heat jump AC{Tq) for 
H=0 is almost absent |25] while AC{Tn) is pronounced. 
However in finite fields of a few Tesla the AFQ HO is 
stabilized and AC{Tq^H) is strongly enhanced. The sta- 
bilization of Aq/ in field is also directly know from RXS 
experiments [8 . This effect will also be present for the 
dynamical resonance. We therefore predict that the res- 
onance peak at R will appear already in the temperature 
range Tn < T < Tq when comparable fields are applied. 
We note that even in the case of a single superconduct- 
ing order parameter the temperature dependence of the 
intensity generally deviates from the BCS MF behaviour. 

In summary the recent INS experiments [1 require a re- 
thinking of the HO phenomena in CeBe. The appearance 
of an itinerant spin exciton resonance at the AFQ wave 
vector proves that the previous restriction to localized 
4f states in CeBe for the hidden AFQ order is oversimpli- 
fied. The neglect of itinerant aspects can no longer be up- 
held. The theory presented here is therefore built on the 
delocalized heavy quasiparticle states. They are gapped 
due to the effect of hybridization and AFQ/AFM type 
particle-hole condensation leading to an enhanced mag- 
netic response at the R-point. Due to quasiparticle inter- 
action a pronounced spin exciton resonance at this wave 
vector appears. Its salient features of momentum, en- 
ergy and temperature dependence are in agreement with 
experimental observation. Therefore CeBe is the first 
non-superconducting heavy fermion example with a spin 
exciton resonance excitation originating in the AFQ hi- 
den order state. 

We thank D. S. Inosov for communicating experimen- 
tal results and M. Siahatgar for useful discussions. 
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